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Abstract 

We consider a DA — D8 — D8 brane construction which gives rise to a large N QCD at 
sufficiently small energies. Using the gravity dual of this system, we study chiral phase 
transition at finite chemical potential and temperature and find a line of first order phase 
transitions in the phase plane. We compute the spectral function and the photon emission 
rate. The trace of the spectral function is monotonic at vanishing chemical potential, but 
develops some interesting features as the value of the chemical potential is increased. 
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1. Introduction 



There are indications that heavy ion collisions at RHIC produce strongly coupled 
quark gluon plasma (sQGP) (see [[1] for a recent review and a list of references). String 
theory gives rise to a number of gauge/gravity dual pairs where strongly coupled gauge 
theory is related to gravity. Hence, one may expect gauge/gravity duality (at finite temper- 
ature) to be particularly useful for studying the qualitative properties of sQGP. Moreover, 
one may hope that gauge theory intuition and direct connection with experiment will bring 
new insights about gravity itself. 

A particularly attractive model where the low energy degrees of freedom are just that 
of large iV c QCD arises from the collection of N c D4 and Nf D8 — D8 branes . (For some 
work in this direction see P-|T3l| ) In this model one of the spatial directions along the D4 
branes is compactified on a circle of radius i?4, which sets the scale of the glueball masses 
and confinement-deconfinement transition in the theory. Another important parameter in 
the theory is asymptotic separation of the D8 and D8 branes, denoted by L. It sets the 
scale of the meson masses and determines the temperature of chiral phase transition in 
the theory |||. In the regime where L <C R4 the two sets of phenomena happen at widely 
separated scales [[|. 

In this paper we restrict the analysis to this regime and extend the results of H 
to nonzero chemical potential u, dual to baryon number. We find the line of first order 
phase transitions in the T — u plane. We also compute the spectral function and the 
photoemission spectrum in the quark-gluon plasma and find an intriguing picture. While at 
small chemical potential the spectral function is featureless, at larger values of a it exhibits 
a characteristic peak at intermediate frequencies. We also observe that the photoemission 
spectrum is redshifted as we increase the value of a. The nature of these phenomena is 
not completely clear to us, and it remains to be seen whether these are generic features or 
they are associated with the particular string mo delS. 

The rest of the paper is organized as follows. In the next Section we review the near- 
horizon geometry of the D4 branes at finite temperature and zero chemical potential. As 
usual, in the limit g s N c ^> 1, N c — > 00, Nf = finite, one can use the DBI action of 
D8 branes to study the physics of fundamental matter. At sufficiently high temperature 



1 The emission of photons associated with the weakly gauged U(1)r symmetry of Af = 4 SYM 
has been previously calculated in [14]. As is well known, in this case nonzero chemical potential 
dual to the U(1)r charge destabilizes the theory. 
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or chemical potential the thermodynamically preferred configuration is the one where D8 
branes fall into the horizon. In this situation baryonic U(l) is no longer broken and the 
corresponding current can be coupled to a photon via the term of the form eJ ^A^ . The 
number of photons emitted per unit volume per unit time is given by 



d 3 k e 



dv= ^wf v °^ mk0 ^ (L1) 



where 



C<{K) = J d A Xe- lKX (J fI (0)MX)) = n b (k°) X ^(K) (1.2) 



In eq. (|1.2|) rif,(k ) = l/(e^ fc ° — 1) and Xfj,u(K) is the spectral function, proportional to the 



imaginary part of the retarded current-current correlation function, 

Xia ,{K) = -2ImC; e J(K) (1.3) 
which we compute using Lorenzian AdS / CFT prescription of []15 . 



In Section 3 we explain how the phase transition takes place at finite chemical poten- 
tial. We generalize the photoemission computation for the case of nonvanishing chemical 
potential in Section 4. Surprisingly, we find that the spectral function exhibits a bump 
whose physical significance is not clear to us. We discuss our results in Section 5. Asymp- 
totic behavior of the spectral function is computed in appendix. 



2. Photoemission at zero chemical potential 

We will consider the near-horizon geometry of the D4-branes in the limit studied in 
0, where the direction transverse to the D8 branes is non-compact. The near-horizon 
geometry of the D4-branes at finite temperature is: 

3 3 

ds 2 = HQ 2 (dxidx* + f(u)dt 2 + (dx^f) + 0Q 2 (J?L + t^dnj) (2.1) 

where f(U) = 1 — U^/U 3 . Here t is the Euclidean time, % = 1, . . . , 3; U and O4 label the 
radial and angular directions of (x 5 , . . . , x 9 ). The parameter R is given by 

R 3 = 7rg s N c l 3 = TVX (2.2) 
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where in the last equality and in the rest of the paper we set a' = 1. The fourbrane 
geometry also has a non-trivial dilaton and RR-four form 

e = 9«[ft) ' F 4 = dC 3 = ~^ e 4 (2-3) 

where Vg and 64 are the volume and the volume form of the unit 4-sphere. Finite temper- 
ature implies that t is a periodic variable, 

t~t + (3 (2.4) 

On the other hand, in order for (|2.1|) to describe a non-singular space, t must satisfy 

4tt# 3/2 

t = t + — 2.5 

Hence, the temperature is related to the minimal value of U, denoted by Ut, as T = 
3U^ /2 /4nR 3 / 2 . 

The induced metric on the D8 branes takes the following form (we assume that the 
temperature is sufficiently high, so the chiral symmetry is restored Q) 

ds 2 = fadsW - f(U)dt 2 ) + QQ " 3/2 ^ + R^U^dQl (2.6) 

In order to find the photon emission rate we need to find current-current correla- 
tors. According to gauge/gravity duality the current- current correlators are computed by 
analyzing the linearized perturbations of the dual gauge field, which in our case is the 
elecromagnetic field living on the D8 brane. In the following we compute photoemission 
due to a single brane. The total result contains an additional factor of 2Nf. 

At non-zero temperature, rotation plus the gauge invariance implies that the correlator 
has the form 

c£* W = ft^W + (Mn L (M (2.7) 

where the non-zero components of the P^ v {k^) are 

P^K) = - ^ (2.8) 
k z 



and 

P*=„ -I 

k 2 
3 



Substituting in (|1.3Q we find that the trace of the spectral function is 



X 



I = -AImIi T {k^) - 2ImU L (k^ 



(2.10) 



In this paper we are interested in computing the rate of production of real photons, hence 
we will only consider lightlike momenta. In this case I1 L (A; M ) vanishes (otherwise the 
correlator ( |2.7|) would be singular on the lightcone) and only Il T (/c M ) will contribute to the 
trace of the spectral function. 

Now we turn on electromagnetic field on the brane. We will be interested in field 
configurations which do not depend on the coordinates on S 4 , we also set Am along the 
S 4 to zero. The DBI action is 



S = -g s T 8 V s R s J d 4 x J dUUe-*y/-det(g AB + 2nF AB ] 



(2.11) 



where (A, B = 0, • • • 3, U), T$ = l/(27r) 8 (? s is the tension of the brane and Vg = 87r 2 /3 is 
the volume of the unit four-sphere. Expanding it in the second order in the field strength 
we find 



S =-(2n 2 )T 8 VsR 3/2 f d 4 x J dUU 5/ ' 



1 (R 



f(U) \U 



E 2 -F 2 u+ 



i<j 



[2.12) 



It is convenient to make a change of variables 

2 ^3/2 



y 



u 



u 



4R 3 



(2.13) 



which brings ( |2.12| ) to the form 



S 



where 



64 



tc 4 T s R 6 / d 4 x / dyy 



f(y) 



if 



i<j 



f(y) = i 

Vt 



2nT 



(2.14) 



(2.15) 



The space is restricted to lie in the region y e (0,i/t), whose upper limit corresponds to 
the black hole horizon. It is convenient to perform Fourier transform in the space-time 

d 4 k , ; 



(2n) 4 
4 



(2.16) 



Choosing k^ 
field 



-u, 0, 0, q), we find the following e.o.m. satisfied by the transverse electric 

.2 _ q 2 f 



-E± = 



(2.17) 



Let us discuss the asymptotics of this equation in the two regimes: y — > yr and y — > 0. In 
the near horizon limit the equation is solvable and E± is 



E ± = (y T -y)±™/ 2 



(2.18) 



where w = utyr/S = u)/2tcT. To compute the retarded correlators we should impose the 
incoming wave boundary condition, which fixes the sign in the exponent of the (|2.18| ) to 
be negative. Now let us discuss the opposite limit y — > 0. In this regime there are again 
two solutions which have the following asymptotics 



Fj = 1 + a 3 y 3 /y 3 , + 
Fn = y 3 /y 3 T + ■■■ 



(2.19) 



Setting as = uniquely determines Fj. The solution of ( |2.17| ) with the incoming wave 
boundary condition has the form 



E± =AF I + BF II , 

where B is proportional to A. 

The boundary action for E± can be determined from 



(2.20) 



SS 



5E ± (y = 0,h^ 



\EOM 



lim 



SS 



128n 4 T 8 R 6 
3u< 



y^o SE' ± (y,k^ 



hmy- 2 f(y)E' ± 



:2.21) 



Using ( |2.2U| ) and varying this equation with respect to E±, we find 

S 2 S , 5 2 S 



B 



5A ± (0,k^)5A ± (0,-k^ 
Hence 



8E ± (0,k,)8E ± (0,-k,) ~ 128 ^ R6 yr 3 A 



4n 



7m(n' J '(^)) = -^{\T)T 2 NJm 

and the trace of the spectral function for lightlike momenta 

16tt 



B_ 
A 



(2.22) 



(2.23) 



X 



£ = -4ImU T (q 



0J 



27 



-XT 3 N r Im 



B 
A 



(2.24) 
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To determine the spectral function we solve the equation ( |2.17|) numerically with the 
incoming wave boundary condition and determine ^ as a function of u (we set q = u). 
The result is depicted on Fig. 1. 

For better numerical convergence near horizon it is convenient to rewrite the equation 
( |2.17| ) in terms of slow varying mode defined as 



E± = (y T - y) • 

Substituting in (gT/j ) we find that (p satisfies 

6y 5 



[2.25) 



<p"(y)+ 



IW 



<p'{y)+ 



1 



6y £ 



2\ iw 1 



W IW * 1 ' '" ' N ' '" 1 ..,.,2. ,6. ,4 

T + T; (y T -y) 2 \y^~y 6 ' yj 2 y T -y ' (yf - y 6 ) 2 



<p{y) = o 

(2.26) 

As expected the exponents of this equation near y = yr are and iw. The solution with 
zero exponent corresponds to E± with incoming wave boundary condition. 




Next we would like to discuss the asymptotic behavior of at high and low temper- 
atures. We start by discussing the low temperature regime ui/t — > <x>. It turns out that 
in this regime we can solve (|2.17|) analytically using generalized WKB approximation. We 
perform this calculation in Appendix A and find 



16tt 



(XT)T 2 N c Im 



B_ 
A 



16n 
~81 



XT^Nr sin — 



37^ ^ F (f) (UVT 



SjT(%)\ 8 



(2.27) 
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Fig 2. 3x^(u;)/(a;j/T) 3 / 4 in the units of ^-\T 3 N C as a function ofuyr- The con- 
stant line represents the approximate solution ( 2.2Tj) . 



This approximation works very well at high frequencies as can be seen from Fig. 2. 

In principle one can improve this result by including higher order terms in the expansion 

in powers of uj~ 1 . 

In the high temperature regime we can solve the equation ( |2.17|) perturbatively. The 
leading behavior is computed appendix A and yields 



16tt 
~2T 



(XT)T 2 NJm 



B_ 
A 



16n 
~81 



\T 3 N c uy T . 



(2.28) 



3. Chiral phase transition at finite chemical potential 

In this section we Wick rotate time coordinate to study the chiral phase transition 
at finite chemical potential. Turning on chemical potential is equivalent to turning on 
non-trivial, purely imaginary A { q\u). Then the average Ad charge density is given by the 
variation of the action w.r.t. Aq°^(oo) 

P = iT — my— \EOM = iT lim — (3.1) 

5A^\U = oo) U ^°°5A'£ 0) (U) 

To find the action we allow non-trivial dependence x 4 = r on the radial coordinate U and 
compute the induced metric on the world volume of the D-brane. The resulting action is! 

S = Tsi y 3/2 J d 3 x j dUU^ 2 \Jl + 4tt 2 F ( ° )2 + (U/R) 3 f(U)r'2 (3.2) 



2 We use A\j = gauge. 
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It is convenient to replace — > iA^ 1 ' , so that fi = A^ 1 ' (oo) + const (we will determine 



,(o) 



(o) ( 



this constant below). Then the equations of motion which follow from ( |3.2| ) are 



d 

dU 

_d_ 

dU 



2nU 5 / 2 A' W 



l-(2 7 r4 0) )2+(5) 3 / (^)r' 2 J 

u ll / 2 f(uy 

l-(27r4 0) )2+(|) 3 /(C/)r' 2 J 



(3.3) 



which should be solved with boundary conditions r(oo) = ±L/2 and Aq (t = ±L/2) = 
u + const. As in the case of zero chemical potential there are two types of solutions: r' = 0, 
which has U(Nf) x U(Nf) chiral symmetry and the U shaped solution, which breaks the 
symmetry down to U(Nf). The phase in which chiral symmetry is broken does not have 
free charges, since the quarks form qq pairs which are neutral (the baryons are infinitely 
heavy in N c —>■ oo limit) . Hence we expect that this phase should be unaffected by non-zero 
chemical potential. On the other hand the phase in which the chiral symmetry is restored 
has free charges and should be sensitive to non-zero chemical potential. We come to the 
same conclusion by analyzing the solutions of ( |3~3| ) . Indeed it is easy to see that for the U 
shaped solution Q3.3|) imply that A^ is a monotonic function of r, hence in order to satisfy 
the boundary conditions one has to have A'q °' ) = oil. The situation is very similar to the 
Hawking-Page transition in AdS space in the presence of chemical potential. In this case 
the thermal AdS solution does not allow non-zero charge density, while the AdS black hole 



does [liU . The reason for this behavior from the CFT side is the same as in our situation: 



the AdS solution corresponds to the confined phase in which we do not have free charges. 
For t' = the action (|3.2| ) takes the following form 



S 



T 8 V S R 3 / 2 



T 



d 3 x / dUU 5/2 \ 1 -4tt 2 F, 



ou 



which yield e.o.m 



27^/2^(0) 



1-4tt 2 (^ 



(0)n 2 



c 



3 The equations ( |3.3| ) were obtained in [11], but their for the U shaped solution is not a 
continuous function at r = 0. 



(3.4) 



(3.5) 
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To fix the constant in the relation between \i and Aq (oo) we require Ad charge density to 
vanish for \x = 0. This implies 



^ = 4 0) (oo)-4 0) (t/ T ) 



(3.6) 



Setting Aq(Ut) = 0, the chemical potential is 



^ = 4 0) (oo 



2tt 



dU 



u T 



C 2 



C 2/5 



C 2 + U 5 2tt 



Ut (I 1 6 _ul 

C 2/5 2 1 I 5' 2 ' 5' C 2 



The thermal expectation value of 4d charge density is related to C as follows 

p = 27iT 8 V s R 3/2 C . 



(3.7) 



(3.8) 



The relation between the chemical potential and the charge density simplifies considerably 



for C 3> U^ 2 and C <C U% ■ For the former we find 



5/2 



V2 m±)r K5 



12tt 2 



6 U "5/2 



7T 



while for the latter 



4n 



N r \~ 1/2 p 5/2 



3 5 p 3 



27 (AT)TV+ 27l37r3 (AT) 



(3.9) 



(3.10) 



Rewriting the conditions of reliability for (|3.9| ) and ( |3.10| ) in terms of /i, A and T, we find 
that (Q is reliable for \i > (AT)T, while ( gig ) for \i < (AT)T. It is instructive to 
compare these relations with corresponding relations for 4d free fermions. In this case the 
charge density is 



Pfre 



Nc_ 
6n 2 



[piv 2 T 2 + p 3 ] 



(3.11) 



We see that while small \i behavior in our system bears some similarity with the free 
fermion case, the large /i behavior is completely different. 

To understand the details of the chiral symmetry restoration we need to compare the 
grand potentials of the two phases 



sn = n st - n v = t 8 v s r 3/2 

) 

dUU 5 ' 2 

U 



U tt5 

dU- 



u T Vu 5 + c 5 

JJ5/2 



+ 



VU 5 + c 5 



1 + 



f(U )U* 



f(U)U8-f(U )U* 



(3.12) 
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where Uq is the value of U at the tip of the U shaped solution. Let us first discuss the 
LT — > case. Solving numerically the equation SQ — we find the critical value for C 

U 



1.83 



(3.13) 



Using (|3.3|) to express Uq in terms of the asymptotical separation L and substituting into 
([T7|) we find 



8R d L~ 2 r(3/10)r(6/5) 



£(9/16) 



0.22L~ A X 



(3.14) 



1.83 v^F 
The phase diagram at non-zero T is given by 

H C {T) = <P(LT)\L- 2 , (3.15) 

where 4>(x) is monotonically decreasing function, which can be determined from ( |3.12| ). 



4. Photoemission at finite chemical potential 

Now we would like to compute the photoemission in the presence of non-zero chemical 
potential. As before we turn on electromagnetic field on the brane and expand fl2.11|) in 
the second order in the field strength. Note that det (gAB + Z^Tab) is an even function of 
T, hence 

det (gAB + ^kFab) = det^s(l + (JF 2 — terms) + (JF 4 — terms)) 



(4.1) 



We set Ao = A { ° ] + A , A B = A B for B ^ 0, T m = F^j + F ou and Tab = F AB for 
AB 0U. In addition to the DBI action there is also contribution coming from the CS 
term 

Scs = m j Tr [exp(27rJ- 2 ) A C 3 ] = ^9^p JtaTATAC 3 (4.2) 
where fig = l/(27r) 8 g s . Integrating this action over the internal S 4 we find 



(0) 



Scs 



'24tt 2 



AAT AT 



(4.3) 



Using the expansion (|4.1| ) we find! 

S = -2tv 2 T 8 V s R 3/2 J d 4 x J 



dU 



Vu 5 + c 2 



U 5 + C 2 fR\ 3 -* 2 (U 5 +C 
i \ E 



f(U) \U 



2\2 



U 5 



F u+ 



i<j i 



+ Scs 



(4.4) 



Note that there is a term in action linear in Fou- This term enforces the equations of motion 

,(0) 

ou 



for F ?) hence it is omitted in the expression for the action. 
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The action in coordinates ( |2.13|) is 



S = - y n 4 T 8 R G I d 4 .r 



dyy~ 



1 + CV7z/t 



i + cV7v* & 



E^-(i + cV°/y 1 T°) 2 Fl+ 



Oy- 



4 + (i + cV M-°)/(!/)£4 



i<3 



+ Scs 



where 



C 2 yf 



C 2 



P 



2 10 R 15 (AT) 4 iV c 2 T 6 
The variation of the Chern-Simons term is 

SS = i^r ! 5 A A A F 
8tv z 



Taking into account that 



dA o 0) -_ 8R s_ c y 2 y~T 



dy 



2nJl + C 2 y 10 /y]9 



we find that the equations of motion for the transverse modes are 



// if 3 5 



uj 2 - 



f y V{1 + C 2 y 10 /y^ 



— / 

(i+cv° / v ™y 



3 H „ yT E 2 = 

(1 + CV°M?)/ 
! _ _ q 2 f 

CjO — 



E"+\ j - + 



f y y (1 + C 2 y w /y^ 



E' + 



3 ~ yT - E x = 



(1 + CV%/T°)/ 



(4.5) 
(4.6) 

(4.7) 



(4.9) 



Note that the asymptotics of this equation as y — > yx and y — > are the same as for ( |2.17|) . 
The solution with boundary conditions ( |2.18|) for E± = E\ and Ei\ y = VT = E' 2 \ y=yT = 
has a form 

E x = A X F! + BxFjj (4.10) 



Ei = AfiFj + B2F11 
11 



(4.11) 



where the coefficients Ai, Bi can be computed numerically. By rotational invariance, 
boundary conditions (|2.18|) for E± = E2 and Ei\ y=yT = E[\y =yT = give rise to the 
solution 

E 2 = AiF I + B 1 F n (4.12) 



E, 



-A2F1 — B2F11 



(4.13) 



To compute x% we need to take the linear combination of the solutions such that the 
leading term in E2 (or E\) near the boundary vanishes, and then repeat (|2.21|) -( |2.23| ). 
The spectral function for lightlike momenta is 



AImIi T {q = uj) 



16n 



(XT)T 2 NJm 



AxB x + A2B2 
A\ + Al 



(4.14) 



Solving (|4.9|) with incoming wave boundary condition numerically we find the spectral 
function in the case of non-zero chemical potential. The result is depicted on Fig 3. 
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Fig 3. 3x£(uj)/(ujy T ) in the units of^XT s N c as a function of uyr for C = 200 



The result is somewhat surprising: while the spectral function for zero chemical potential 
is featureless, we see that at non-zero chemical potential things change and we get a second 
maximum. The nature of this peak is not clear to us. 

The leading behavior of the spectral function at high temperature is computed in 
appendix A 

^ ^(\T)T 2 N c uy T ^l + C* . (4.15) 



X 



16tt 
~2T 



XT 3 NJm 



A1B1 + A 2 B 2 



12 



5. Discussion 



In this paper we extended the results of [|J to the case of nonvanishing chemical 
potential. We found a line of first order phase transitions in the \x — T plane which 
intersects the axes at (fx = 0,T = 0.15L -1 ) and (/i c = 0.22L _2 A,T = 0). Note that 
the value of \i c is of the order of constituent quark mass in the phase with broken chiral 
symmetry!. 

We also studied photoemission at finite temperature and chemical potential. At small 
chemical potential the spectral function is featureless while at larger values of /i it exhibits 
a characteristic peak at intermediate frequencies. 

Using (|1 . 1|) we rewrite the photon emission rate as function of the frequency 



dT _ e 2 ux» 



du 47T e w / T - 1 



(5.1) 



Then for the zero chemical potential we plot the spectrum on Fig 4. using the results of 
section 2. For comparison we also plot the black body radiation spectrum on the same 
graph 

dT P u 2 



" * (5-2) 



dw &>I T - 1 




2 4 6 8 



Fig 4. dT/duo in the units of ^-{\T)T*N C as a function of uj/T-red curve and the 
black body radiation spectrum (in arbitrary units)-blue curve. 



The constituent quark mass is given by the mass of the string which goes from the tip of 
D8 brane all the way down to the horizon. In this system it is much larger then the meson mass. 
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Fig 5. 



Vi+c 2 



dr 



in the units of -^-(AT)T N c as a function of u/T. Red, green and 



blue curves correspond to C = 1000, 100, 10 respectively. 



We observe that at zero chemical potential the spectrum of radiation is very similar to the 
Planckian spectrum. As we turn on the chemical potential things start to change. The 
spectrum becomes more and more red-shifted as we increase the chemical potential. 

The natural question is whether something similar happens in other models. One 
particularly simple way of incorporating fundamental matter involves adding D7 branes 



to the stack of D3 branes [ PT|j . The resulting theory is four dimensional, but the presence 
of bosonic fields charged under the U(l)v does not allow turning on nonzero chemical 
potential!. Similar problems persist for the D4-D6 system recently studied in |18|jB. 

The model we consider does not have charged bosonic fields. On the other hand 
one must bear in mind that it involves four dimensional fundamental matter and five 
dimensional gauge field and only becomes four dimensional QCD for energies much smaller 
then l/i?4. This is the reason for the appearance of the factors (AT) in the expressions for 
the spectral functions in this paper. (In comparison, such factors are absent in the four- 
dimensional model of ||14|| .) It would be nice to better understand the physical origin of the 
features of the photoemission spectrum in strongly coupled gauge theories at nonvanishing 
chemical potential. 

Acknowledgements: 



6 One may still induce nonzero charge density by condensing the corresponding bosonic fields. 

7 The value of u in this system is bounded from above by the scalar mass, while we observed 
that the nontrivial features in the spectral density appear at very large values of /i. 
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Appendix A. Asymptotics for the spectral function 



In this appendix we find the asymptotic behavior of the spectral function at high 
and low temperatures. We start by discussing the low temperature regime ivyr — > oo. 
It turns out that in this regime we can solve ( |2.17|) analytically using generalized WKB 
approximation. We will follow the approach of ||20||. First we introduce new variable 



E\ 



-E 



(A.l) 



to eliminate the first derivative term in ( |2.17|) i 

2y 12 - 13y 6 + 2 



E'l + 



uj 2 y 6 



E ± = 



(A.2) 



(l- y 6)2 y2 (1 _ 2/ 6)2 

In the uj — > oo regime the first term in the brackets dominates almost everywhere on the 
segment y G (0, 1). As we approach y = point the second term starts contributing. Keep- 
ing only the most singular contribution from the second term we obtain the approximate 
equation 



E'i + 



u 2 y 6 



Rewriting 



;i-y 6 ) 2 

in terms of new variables 



dy 

dz 



V s 



E ± =0 



(A.3) 



1-1/ 



6^2 



w 



(A.4) 



we obtain 



d 2 W 
dz 2 



+ 



u 2 z 6 - 



W = 0. 



This equation can be solved analytically in terms of Bessel functions 



3 

(00\ 8 



W(z) = AT [~ ^ y/zJ_z{uz*/4) + BY — ^ ^z~Js (uz 4 /4) 



11\ /u 



(A.5) 



(A.6) 



Here and below we rescaled y — > y/yr and uj —> ojyr for simplicity. 
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As we will see shortly the first term corresponds to the Fj in ( j2.2Q|) , while the second term 
corresponds to Fjj. Hence we identified the integration constants with A and B. Then 
the approximate solution to (|2.17|) is 



To find the spectral function we will need the behavior of z as y — > 1 and y 
( |A.4| ) we find 



(A.7) 
0. Using 



z 4 /4 



1 



2^3 V 



arctan 



l + 2y 

v/3 



— arctan 



2y-l 

J 



- ^log(l-y)- 

6 



^ log(l + y) + — log(l - y + y 2 ) + — log(l + y + y 2 ) 



TV 



(A.t 



Hence 



z =y + 



for y — >• 



z = ~ g log(l - y) H , for y 

Using the asymptotic behavior of the Bessel functions 



(A.9) 



1 



X 



for x —>■ 



(A.10) 



r(a + 1) \2. 

we confirm the identifications made in ( |A.6| ). To find the spectral function we should look 
at the solution ( |A.7| ) in the opposite regime y — > (z —>■ oo) and impose the incoming 
wave boundary conditions 



E 



± 



(!-»)■ 



;/6 



g 

~ e 4 



Using the asymptotic form of Bessel functions at large value of the argument 

Y 



J a (x) 




COS X 



we find 



Ira 



B 
A 



TlX 



sin 



an tt 



3tt 



r(f) ( 



3 



(A.11) 



(A.12) 



(A.13) 



sjr^i) Vs. 

The case of the finite chemical potential can be treated similarly, although the calculation 
is more involved and we will not pursue it here. 

Now let us discuss the high temperature case to — ► 0. In this regime we can solve 
equations ( |2.17|) and the corresponding equation for the finite chemical potential 
perturbatively in u = 3w 



E ± (y) = (1 - y)-^V(y) = (1 - y)-^(F^(y) + wF^(y) + ■ 



(A.14) 



16 



We start by discussing the zero chemical potential case. Then <p(y) satisfies equation ( |2.26| ). 
In the w° order the equation has two solutions one of which is constant while the other 
is logarithmically divergent at the horizon. The solution with incoming wave boundary 
conditions correspond to the constant solution which we set to one. Substituting (|A.14| ) 
into ( |2.26| ) we find that F^(y) at w order satisfies 



6y 



2 

+ - 

y 



6y 



5 2 
+ 



(A.15) 



.(! - y) 2 V 1 - y 6 y/ 1 - y. 

Solving ( |A.15|) we find 

F« = ^ ((5 - Md) log(l - y) - (1 + 2iC t ) log ( ^^l^f X) (A. 16) 

The coefficient in front of the first term should be zero to satisfy the boundary condition. 
Using this explicit solution we find that near horizon E±_ is 



E±(y) = 1 + iwy 3 + 



(A.17) 



Hence to the leading order in w 



Im 



B_ 
A 



w. 



In the case of finite chemical potential we have 



E 1 = (1 - y)-r<pi{y) = (1 - y)-T{Fi°\y) + wF^( y ) + ■ ■ 
E 2 = (1 - y)-* V2 (y) = (1 - y)-^(F!f\y) + wF^\y) + • • 



(A.18) 



(A.19) 



The incoming wave boundary conditions for Ei and E 2 \y= yT = E 2 \ y=yT = imply that 



Fi 0) (y) = l and F 2 ^> (y) = 



(0), 



(A.20) 



Then in the next to leading order the two equations in (|4.9|) decouple and we find that 
F^ satisfies the following equation 

6y 5 



3 5 
- + -- 



i - y b y yi + c 2 y 



10 



Fi 1} (y)'+ 



l 



6y E 



+ 



i-y 6 y yi + c 2 y 10 J i-y 



(A.21) 
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satisfies similar equation, which we will not need at this order. Solving QA.21 ) we find 



Fi 1 \y) = C 1 f- 
Jo (1 



-Ay + ilog(l - y) 



(A.22) 



{l-y*)dl + C*y 



,10 



The integral above is logarithmically divergent near y = 1. Imposing the regularity condi- 
tion we find C\. Near y = we can again approximate the integral by dropping the square 
root in the denominator. Hence we find 



Ex{y) = l + iw\/l + C 2 y 3 + 
E 2 (y)=wF^\y) + 



(A.23) 



and 



Ira 



A 1 B 1 +A 2 B 2 



Im 



A X 



W 



1 + C 2 . 



(A.24) 
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